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Abstract 



The nature of the transition from quantum tunneling at low temperatures to ther- 
mal hopping at high temperatures is investigated in a scalar field theory with cubic sym- 



metry breaking. The bounce solution which interpolates between the zero-temperature 



and high-temperature solutions is obtained numerically, using a multigrid method. It 
is found that, for a small value of the symmetry-breaking coupling /, the transition is 



first-order. For higher values of /, the transition continues to be first-order, though 
weakly so. 

I. INTRODUCTION 



On 
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The problem of decay of a metastable state via quantum tunneling has important applications 
in many branches of physics, from condensed matter to particle physics and cosmology. In 
the semi-classical approximation, the decay rate per unit volume is given by an expression 
of the form 

r = A e- SE , (1) 

where Se is the Euclidean action of the bounce: the classical solution of the equation of 
motion with appropriate boundary conditions. The bounce has turning points at the con- 
figurations at which the system enters and exits the potential barrier, and the analytic 
continuation to Lorentzian time at the exit point gives us the configuration of the system at 
that point and its subsequent evolution. The solution of the equation of motion looks like 
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a bubble in four dimensional Euclidean space with radius R and thickness proportional to 
the parameter in the symmetry breaking term. If there is more than one solution satisfying 
the boundary conditions, the one with the lowest Se will dominate Eq. ([!]). The prefactor A 
comes from Gaussian functional integration over small fluctuations around the bounce. The 
formalism at zero temperature is well known Jl|, [|; the least action is given by the bounce 
which is 0(4) invariant [H]. 

At nonzero temperature, the above formalism has to be modified, since the metastable 
state can decay due to classical thermal motion. In the context of quantum mechanics, the 
extension of the formalism to finite temperature was carried out by Affleck . He argued 
that a second order phase transition from the quantum to the thermal regime takes place at 
some critical temperature. Chudnovsky j5| showed that this is not always true. Depending on 
the shape of the potential barrier, the crossover from thermally assisted quantum tunneling 
to thermal activation at high temperature is either first-order (i.e., the first derivative of T 
at the transition temperature is discontinuous) or second-order (i.e., the first derivative of T 
is continuous, but the second derivative is discontinuous at the transition temperature). 

In the context of quantum field theory, little work has been done to explore the nature 
of the transition from the quantum regime (zero-temperature) to the thermal regime (high- 
temperature). Linde suggested that periodic bounces could smoothly interpolate between 
these two regimes. In this picture, at zero temperature the solution is an 0(4) symmetric 
bubble with a radius R. Up to T ~ (2i?) _1 , we have periodic, widely separated bounces, 
and beyond this temperature they start merging into one another producing what is known 
as a "wiggly cylinder" solution. As one keeps increasing the temperature these wiggles 
smoothly straighten out and the solution goes into an 0(3) invariant cylinder (independent 
of Euclidean time r), which dominates the thermal activation regime. 

Garriga |7| extended the work of Chudnovsky by studying bubbles in the thin-wall ap- 
proximation (TWA) at high and low temperatures, as well as the wiggly cylinder solution 
at intermediate temperatures. He showed that in the TWA the transition is always first- 
order. However, although motivated by quantum field theory, this work wan not truly 
field-theoretical. 

Ferrera || carried the investigation beyond the domain of validity of the TWA. In fact, 
this was the first fully field-theoretic investigation of bubble formation at arbitrary tempera- 
tures. (Such a study must solve the equation of motion at intermediate temperatures, which 
requires solving a partial differential equation.) Ferrera studied a model field theory which 
was if 4 theory with a symmetry-breaking term of the form ftp. His result is that only for 
very large wall thickness (i.e., large /) a second-order phase transition takes place, while for 
all other first-order phase transition occurs. 
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In an earlier paper |J, we studied phase transitions in two kinds of p 4 theory, with 
symmetry-breaking terms proportional to (p and p 3 respectively. We obtained accurate 
numerical solutions to the equation of motion in the zero- and high-temperature limits, and 
found that, even for a fairly large value of the p 3 coupling, the thin- wall approximation 
(TWA) holds to a few percent. An analytical solution for the bounce was obtained, which 
reproduces the action in the thin-wall as well as the thick-wall limits. We also investigated 
the dependence of T*, defined by S4/S3 = l/T*, on the dimensionless symmetry breaking 
parameter /. For small /, the TWA behavior T* oc / is obtained, while at larger values of 
/ there is a smooth departure from this || . 

It has been argued recently that the leading correction to the tree potential due to one- 
loop effects at finite temperature is proportional to p 3 rather than tp []10[ . It therefore 
becomes natural to investigate p 4 theory with p 3 symmetry-breaking at arbitrary tempera- 
tures. In the present paper, we undertake such an investigation. Our aim is to obtain the 
bounce solutions at all temperatures, to compute the action and to determine the nature 
of the phase transition from quantum tunneling to thermal hopping. Our work therefore 
follows a parallel track to that of Ferrera || , and it is of interest to compare the results for 
the two theories. We find that, for small values of /, the action shows a kink when plotted 
against temperature or its inverse f3. This can be identified as a first-order transition. For 
larger values of / the transition is weakly first-order. The second-order behavior seen by 
Ferrera for / = 0.75 is not seen here. 

In Sec. II we review the formalism for bounce solutions at finite temperatures. In Sec. Ill 
we discuss our numerical algorithm. In Sec. IV we present our results. Finally, conclusions 
are presented in Sec. V. 

II. FINITE TEMPERATURE BOUNCE 
SOLUTIONS IN FIELD THEORY 

Let us consider a scalar field theory with a Lagrangian density 

C{p) = \{d»p) 2 -U{p), (2) 

where the potential U(p) has two minima at p- (false vacuum) and p + (true vacuum). 

In the semi-classical approximation the barrier tunneling leads to the appearance of 
bubbles of a new phase with <p = p + . To calculate the probability of such a process in 
quantum field theory at zero temperature, one should first solve the Euclidean equation of 
motion : 

aA ^«JM, (3) 
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with the boundary condition <p — > as x + r — ► oo , where r is the imaginary time. The 
probability of tunneling per unit time per unit volume will be given by 



.4 



-Se[<p] 



(4) 



where Se[<p] is the Euclidean action corresponding to the solution of Eq. (|) and given by 
the following expression : 



Se[<p] 



d 4 x 



(5) 



It is sufficient in most cases to restrict ourselves to the 0(4) symmetric solution ^(a^+r 2 ), 
since it is this solution that provides the minimum S4 of the action S^fy?] |§. In this case 
Eq. (^|) takes the simpler form 



d 2 if ^3 dip 
dp 2 p dp 

where p = ^/x 2 + r 2 , with boundary conditions 



dU{ip) 
dip 



if — > <y?_ as p 



00, 



dip 
dp 



at p = 0. 



(6) 



(7) 



Now, let us consider the finite temperature case. Following ||, in order to extend the 
above-mentioned results to nonzero temperature, T 7^ 0, it is sufficient to remember that 
quantum statistics of bosons (fermions) at T 7^ is equivalent to quantum field theory 
in the Euclidean space-time, periodic (anti-periodic) in the "time" direction with period 
(3 = T _1 . One should use the T-dependent effective potential U(<p,T) instead of the zero- 
temperature one U((p) = U(<f, 0). Instead of looking for 0(4)-symmetric solution of Eq. (|3|), 
one should look for 0(3)-symmetric (with respect to spatial coordinates) solutions, periodic 
in the "time" direction with period (3 = T~ l . At sufficiently large temperature compared to 
the inverse of the bubble radius R at T = 0, the solution is a cylinder whose spatial cross 
section is the 0(3)-symmetric bubble of new radius R{T). In this case, in the calculation 
of the action Se[^p], the integration over r is reduced simply to multiplication by T _1 , i.e., 
SeW\ = T^ 1 £3 [</?], where S 3 [ip] is a three-dimensional action corresponding to the 0(3)- 
symmetric bubble: 



d A r 



-(Vy?) 2 + f%,T) 



To calculate S%(ip) it is necessary to solve the equation 

dU(<p,T) 



d 2 ip 2 dip 
dr 2 r dr 



dip 



with boundary conditions 



as r 



00, 



dip 
dr 



at 



0. 



(9) 



(10) 



where now r = The complete expression for the probability of tunneling per unit time 
per unit volume in the high-temperature limit (T >> is obtained in analogy to the one 
used in || and is given by: 

T(T) = AiT) e ~ s ^ T]/T . (11) 

At intermediate temperature the action will be given, under the assumption of 0(3) 
symmetry, by 



S E [<p,T\ = 4tt J dr J dr r 5 



and the equation of motion will be 



d 2 cp 2dcp _dU((p,T) 



dr 2 dr 2 r dr 



dip 



with boundary conditions 



(12) 



(13) 



(p — >• <f- as r — ► oo, dip/dr = at t = ±{3/2,0 



(14) 



where (3 is the period of the solution. 



III. ALGORITHM 



The aim of this paper is to investigate the nature of the phase transition for tp 4 theory with a 
ip 3 symmetry-breaking term and compare our results with those of Ferrera M who discussed 
a ip symmetry-breaking term. 

We start with the following Euclidean action inspired by recent work on temperature- 
dependent corrections to the tree level potential [HI : 



Sf = 4 7r dr dr r 



l(^) 2 +(^r±u(ip) 



(15) 



where U(ip) is 



. We rescale the Lagrangian density as 



(16) 



ip^ip/li, r^r^\{i 2 , f = 



(17) 



to get 



Se = -T- / dr / s 'li' i' 2 



l(^) 2 + (^) 2 + l(v 2 -i) 2 -fV 
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The only free parameter in the Lagrangian is /. Taking different values of / will give us 
different shapes of the potential. Hence we will be able to explore the nature of the phase 
transition for thin as well as for thick walls. The equation of motion now becomes 

j£ + f£ + H£ = 2fc»- 9 )-sM (19) 
ot z or A r or 

with boundary conditions 

y?^<^_ as r^oo, d^pjdr = at t = ±{3/2,0. (20) 
To solve Eq. ([19|) we have used a multigrid algorithm. Multigrid methods were introduced 



in the 1970s by Brandt More details are available elsewhere 1121 , |13 |. 

We have used full weighting injection for the projection operator and second-order poly- 
nomial interpolation. For the thin wall we have used fourth order polynomial interpolation 
since more points are coupled to each other than in the thick-wall case. Since we used the 



multigrid algorithm (not the full multigrid [121), we have to start from the finest grid and go 



down to coarse grids and then up and down till the solution converges. Our finest grid has 
513x513 points. We have restricted ourselves to three grids (513x513, 257x257, 133x133). 

As we can see Eq. ([B3|) has a singularity at r = 0. One way to handle it by using an 
arbitrary regularization 1/r — > l/(r + e) with e chosen suitably small. This is the method 
used by Ferrera | 14j . We use an alternate method based on L'Hospital's rule, replacing 



(l/r)d(p/dr at r = by d 2 (p/dr 2 . 

We have started with the zero temperature solution as an initial guess solution, as the 
zero-temperature equation of motion is one- dimensional and straightforward to solve . We 
used this one-dimensional solution after putting it in two-dimensional form. After computing 
the zero-temperature solution by the multigrid method we used it as a guess solution for 
our first finite-temperature solution. To compute the next finite-temperature solution we 
used the previous one as input, and so on. We followed this procedure (i.e., continuation in 
temperature) to cover the whole range of temperatures. 

At zero temperature we have used point Jacobi relaxation and at finite temperature we 
have used line Newton Jacobi relaxation in the r direction because as the temperature is 
increased we have to reduce the step size in the r direction and this means that more points 
in the r direction will be coupled than in the other direction. 

To solve the equation of motion, we need the value ipx c of the bounce at the center. 
For this, we chose a trial value tp c and kept it fixed through the relaxation process (since 
otherwise the process is bound to diverge ||). In the zero- and high-temperature limits, 
ipT c was fixed very accurately by solving the one-dimensional equation of motion using a 
predictor-corrector method and shooting in (p c |§. For multigrid, we have found the same 
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behavior as Ferrera i.e. if we give too high value of ip c the multigrid solution diverges 
while for low enough values it converges. If ip c — ifxc ~ A, the solution converges to an 
overall accuracy level of order A. We found A ~ 0.1% which corresponds to a low level of 
error for the solution. The shooting in (p c , together with the variable amount of relaxation 
required, makes the program very CPU-intensive. A typical running time for a single value 
of temperature is about an hour and a half on an 80586 Pentium machine running at 133 
MHz. 

A drawback of the multigrid method is that we cannot find the solution at an arbitrary 
value of temperature. We have a solution only at discrete values of (3. A constant step 
size in (3 implies the interval between successive temperatures increases as we go to higher 
temperatures. On the other hand, the (3 step cannot be made very small because of the 
reasons stated above. 

Finally, to compute the action for each solution we have used a two-dimensional Simpson 



method that can be found in 15 



IV. RESULTS 

We have used the algorithm described in the above section for solving the equation of motion 

<9V d 2 (fi 2 Old , , , „ 9 . , 

^ + ^ + -^ = 2(^-„)-3V. (21) 

To be able to compare our results with Ferrera ||, we have chosen the same three values of 
the parameter /. The first value is / = 0.25. From [§],/ = 0.25 is within the domain of 
validity of the TWA Fig. 1 represents the shape of the potential for / = 0.25. If we compare 
the action of the zero- and high-temperature solutions with the one-dimensional one, we 
find a difference of order 0.2%. At intermediate temperatures the highest difference is about 
0.9%. Fig. 2 shows a plot of S(T) vs (3. From it, the value of (3 at which the transition takes 
place is approximately 21.15 in our dimensionless units. Garriga has shown that, in the 
TWA, the period at which the transition takes place is given by 

(3c = ^R, (22) 

where R is the radius of the one-dimensional zero-temperature solution. We have measured 
the radius from the center of the bounce to the middle of the wall where the absolute first 
derivative of the field is maximum. Using Eq. ([22|), we obtain a value of (3 C ~ 21.205. The 
difference between the numerical result and the analytic one is about 0.2%, which lies within 
our estimation of errors. Clearly it is a first-order phase transition and this will always 
happen as long as the thickness of the wall is much smaller than the radius of the nucleated 
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bubble. Fig. 3(a) shows the zero-temperature solution. We can see that the wall thickness 
is much smaller than the radius. Fig. 3(b) shows the solution at (3 = 25.3. We can see that 
the two bubbles are far from touching each other. Fig. 3(c) shows the solution at f3 — (3 C . 
In this case we can see that the two bubbles just fit within the period. Finally, Fig. 3(d) 
shows the solution at f3 < (3 C . 

The second value is at / = 0.55. Fig. 4 shows the shape of the potential. It is clear from 
the figure that the TWA is not valid since the energy difference between the false and true 
vacua is larger than the hump at the origin. Fig. 5 shows a plot of S(T) vs (3. From the 
figure, we can see that the action is constant till j3 = 14.4 and after that it starts decreasing 
slowly till it matches with the (3 S3 curve, where S3 is the high temperature action. The 
transition point is at f3 — 9.3. If / = 0.55 were within the TWA, then the transition point 
would be at /3* = 9.9. We can see that the transition point is at a period less than Fig. 
6(a) shows the zero-temperature solution. Fig. 6(b) shows the solution at (3 — 14.4, where 
the two bubbles have not started touching each other. Hence we have constant action till 
that period. Fig. 6(c) shows the solution at f3 — 9.3. We can see that the solution looks 
like a wiggly cylinder but it is far from a static one, so we still do have a first-order phase 
transition, but more weakly than for / = 0.25. Finally, Fig. 6(e) shows the solution at 
(3 < 9.3. 

Finally, Fig. 7 shows the shape of the potential for / = 0.75. Fig. 8 shows a plot of 
S{T) vs (3. From the figure we can see that the action is constant till (3 = 10.3. After that it 
starts decreasing till it crosses the (3S 3 curve at f3 — 6.25. Again, if / = 0.75 were within the 
TWA, then the crossing point would be at f3* = 7.45 . The transition point is at f3 — 6.25, 
which is again less than f3+. Fig. 9(a) shows the zero-temperature solution. Fig. 9(b) shows 
the solution at (3 — 10.3; we can see that wiggly cylinder solution has not started forming 
yet. Fig 9(c) shows the solution at the transition point (3 = 6.25. We can see that the 
wiggly cylinders are formed and still far from a static solution. Hence we have a first-order 
phase transition but more weakly than for / = 0.55. Finally, Fig. 9(d) shows the solution 
at (3 < 6.25. 

V. CONCLUSIONS 

In this paper we have tried to find out whether the phase transition from the quantum 
regime to the high-temperature regime is of first or second order for y? 4 theory with a y? 3 
symmetry-breaking term. We have investigated three values of /. We find that for / = 0.25 
we have a first order phase transition. For / = 0.55, the transition is still first-order but 
weaker than for / = 0.25. Another feature is that the transition point (3 C is less than the 
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point /?* where the zero-temperature solution will match the high temperature one. The final 
value is f — 0.75. Here we have a weak first order phase transition. Again f} c < f3±. In this 
case the departure is more than for / = 0.55. So, up to / = 0.75 we have a first-order phase 
transition, which is different from what Ferrera || found for / = 0.75. We conclude that, 
for the potential with a <p 3 symmetry-breaking term, a true second-order phase transition 
is not seen even at / = 0.75. This can be understood as follows. A second-order phase 
transition only occurs when the hump in the potential disappears. In [[| we found that for ip 
symmetry breaking the hump in U ((f) disappears for / > 0.77, while for if 3 it disappears only 
asymptotically in /. The two forms of the potential can be mathematically transformed into 
each other ||16|| . However, this obscures the physical meaning of the asymmetry coefficient 
/, so we work with the untransformed potential. 

We find that the departure of the transition point (f3 c ) from the point where 0(4) and 
0(3) match is higher for higher values of /. This is expected since, in the domain of 
validity of the TWA, the surface tension of the bubble is independent of the temperature 0; 
hence we must have (3 C = (3*. Beyond the TWA, on the other hand, we can expect departure 
from this equality. 

Finally, in this paper we have studied a zero-temperature potential motivated by the 
finite-temperature effective potential fID |. We have incorporated the temperature only in the 
boundary conditions. In an exact calculation, of course, we should use the full temperature- 
dependent effective potential. However, we expect that there will be a range of temperatures 
and parameter values for which the potential of Eq. (|i~6|) will be a good approximation to 
the full potential. Hence, our conclusions about the nature of the transition to the thermal 
hopping may have relevance to actual physical situations, in particular inflationary cosmology 
and baryogenesis. 
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Figure Caption 



FIG. 1. Shape of the potential for / = 0.25. 

FIG. 2. Temperature dependence of the Euclidean action: S(T)vsf3 for / = 0.25. 
FIG. 3. Bounce solution for / = 0.25. (a) T = 0, (b) (3 = 25.3, 
(c) P = 21.15, and (d) high temperature. 
FIG. 4. Shape of the potential for / = 0.55. 

FIG. 5. Temperature dependence of the Euclidean action: S(T)vs{3 for / = 0.55. 
FIG. 6. Bounce solution for / = 0.55. (a) T = 0, (b) (3 = 14.4, 
(c) (3 = 9.3, and (d) high temperature. 
FIG. 7. Shape of the potential for / = 0.75. 

FIG. 8. Temperature dependence of the Euclidean action: S(T)vs{3 for / = 0.75. 
FIG. 9. Bounce solution for / = 0.75. (a) T = 0, (6) (3 = 10.3, 
(c) (3 = 6.25, and (d) high temperature. 
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U(q>) for/=0.25 
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Figure 2: 
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(a) f=0.25, Temperature= 




Figure 3: 
Figure: 3 (a) 
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(b) f=0.25, Period= 25.3 




Figure: 3 (b) 
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(c) f=0.25, Period= 21.15 




Figure: 3 (c) 
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(d) f=0.25, Thermal hopping 




Figure: 3 (d) 
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Figure 4: 



18 




19 



(a) f=0.55, Temperature= 




Figure 6: 
Figure: 6 (a) 
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(b) f=0.55, Peroid= 14.4 




Figure: 6 (b) 
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(c) f=0.55, Peroid= 9.3 




Figure: 6 (c) 
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(d) f=0.55, Thermal hopping 




Figure: 6 (d) 
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U(q>) for/=0.75 
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Figure 7: 
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(a) f=0.75, Temperature= 




Figure 9: 
Figure: 9 (a) 
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(b) f=0.75, Period= 10.3 




Figure: 9 (b) 
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Figure: 9 (c) 
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(d) f=0.75, Thermal hopping 




Figure: 9 (d) 
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